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Abstract 

"A moving body will come to rest as soon as the force pushing it no longer acts on it in the manner 
necessary for its propulsion. " 

(Aristotle) 



We discuss the transport of a tracer particle through the Bose Einstein condensate of a Bose 
gas. The particle interacts with the atoms in the Bose gas through two-body interactions. In the 
limiting regime where the particle is very heavy and the Bose gas is very dense, but very weakly 
interacting ("mean- field limit"), the dynamics of this system corresponds to classical Hamiltonian 
dynamics. We show that, in this limit, the particle is decelerated by emission of gapless modes into 
the condensate (Cerenkov radiation). For an ideal gas, the particle eventually comes to rest. In an 
interacting Bose gas, the particle is decelerated until its speed equals the propagation speed of the 
Goldstone modes of the condensate. This is a model of "Hamiltonian friction". It is also of interest 
in connection with the phenomenon of "decoherence" in quantum mechanics. It is based on work 
we have carried out in collaboration with D Egli, IM Sigal and A Soffer. 
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I. INTRODUCTION 



We review some recent results on the effective dynamics of particles moving through a dispersive 
wave medium, such as the Bose-Einstein condensate of a very dense quantum gas at very low 
temperature. The motion of the particles exhibits friction as long as their speed is larger than 
the propagation speed of waves in the medium. Friction is caused by the emission of Cerenkov 
radiation. This phenomenon provides an example of recent mathematical results in transport theory. 
A key problem in transport theory consists in deriving the effective dynamics of experimentally 
observable degrees of freedom interacting with unobserved degrees of freedom of a macroscopically 
large environment from an underlying fundamental (Hamiltonian or unitary quantum) dynamics of 
the entire system; as indicated in the table below. 




Fundamental Dynamics 


Effective Dynamics 


Hamiltonian dynamics of point 
particles and waves - celestial 
mechanics; 

Schrodinger equation, mechanics of 
atom gases;... 


Dissipative dynamics of particles 
moving through a medium; 

friction, diffusion, Brownian motion, 
(quantum) Boltzmann equation, dy- 
namics of viscous fluids,... 



The key problem alluded to, above, concerns the passage from the left to the right column in this 
table. Typically, one might study a system consisting of an observable tracer particle, or a whole 
tribe of such particles, coupled to unobserved or only crudely observed degrees of freedom of a 
medium such as the electromagnetic field in an optically dense material, sound waves in a crystal or 
the Goldstone modes of an atom gas exhibiting Bose-Einstein condensation. At low temperatures, 
the emission of Cerenkov radiation by the particle(s) into the medium is the basic mechanism 
causing friction of the particles' motion. When the medium is, initially, in thermal equilibrium at 
some temperature T > 0, it exhibits thermal noise that causes a diffusive type of motion of the 
particle(s) around its (their) mean trajectory. 

The problem of making these remarks precise by identifying the correct effective dynamics of the 
particle(s) and analyzing its properties leads one to study a variety of interesting, novel phenomena 
in the realm of non-linear Hamiltonian evolution equations and quantum dynamics for systems with 
infinitely many degrees of freedom. 



A Dynamics of a particle moving through a medium 
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II. SOME SIMPLE MODEL SYSTEMS 
A. Dynamics of a particle moving through a medium 

We consider a classical particle moving in physical space E 3 . The Hamilton function of this system 
is given by 

H(X,P) = ^ + V(X) (1) 

with IgE 3 the position of the particle, P G R 3 its momentum, M its mass, and V the potential 
of an external force. An example of a potential is 

V{X) =-F ■ X (2) 

where F G M 3 is a constant external force pushing the particle; e.g., the gravitational force close to 
the surface of the earth. 

The particle is assumed to interact with a wave medium. We consider the example of a Bose 
gas exhibiting Bose-Einstein condensation. We would like to treat the entire system within the 
framework of classical Hamiltonian dynamics. This forces us to pass to a limiting regime of the 
system, the so-called mean field limit: The average density of the Bose gas is given by pn '■= Np, 
p > 0, the strength of two-body interactions among the atoms of the Bose gas is Ajv := N _1 X, 
A > 0; we consider the limiting regime where N — > oo, with p and A kept fixed. In this limit, 
the state of the gas can be described by a Ginzburg-Landau order parameter, ip(x) G C, where 
N\tp(x)\ 2 is the density of the gas at the point x G E 3 . In microscopic units, the mass of a 
quantum-mechanical tracer particle moving through the Bose gas is given by NM (i.e., it becomes 
very heavy, as N — > oo) and the external potential is given by NV(X), with M and V fixed. The 
limit N — > oo corresponds to a "classical limit" in which the quantum dynamics of the system 
approaches some classical Hamiltonian dynamics. 

The Ginzburg-Landau order parameter of the Bose gas, ip(x), and its complex conjugate, ip{x), are 
interpreted as complex coordinates of an infinite-dimensional affine phase space, Tbg '■= H^M 3 ), (a 
certain complex Sobolev space). The standard symplectic structure on Tbg leads to the following 
Poisson brackets : 

WWMv)} = {$(x)Mv)} = o (3) 

{V>(x),^(y)} = iS(x-y) 

Quantum mechanically, ip and ip would be interpreted as annihilation- and creation operators, and 
the Poisson brackets are replaced by commutators, [•, ■], with [-, •] ~ — •}. 

In the mean-field limit, N — > oo, the following Hamilton functional is appropriate to describe a 
non-relativistic Bose gas: 

Hbg(A$) = J d 3 x ^L|VV(x)| 2 + a| d 3 ymx)\ 2 -p)Hx-y)(\iP(y)\ 2 -p)] (4) 



B A list of models of decreasing difficulty 
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where $ is a two-body potential of short range and, for reasons of thermodynamic stability of 
positive type, A > is a coupling constant and p > is an average (rescaled) density. We consider 
a regime where the [/(l)-symmetry 

ip(x) i—)- e i9 ip(x), 4){x) e~ ie %b{x), 6£R (5) 

is spontaneously broken. It is convenient to impose symmetry-breaking boundary conditions 

% (x) = ^p+ { p\x) (6) 

with 

P(x) — > 0, as \x\ — > oo 



The Poisson brackets between f3 and (5 are the same as those in ([3]). In terms of the variables f3 
and /3, the Hamilton functional is given by 
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n BG (0,0) = / d 6 x —\vp( x )\*+ 

J L ZVTl ^rj^ 

A J d 3 y (\/3(x)\ 2 + 2^p WW) *(* ~ V) (l/%)| 2 + 2VP StP{y))\ 

The interaction between the tracer particle and the Bose gas is given by 
Hl(X; 0, 0) :=g f d 3 x W(X - x) (\^{x)\ 2 - p) 



(8) 

=g / d 3 x W(X - x) (|/3(x)| 2 + 2^p UP{x)) 



where X is the position of the tracer particle, W is a spherically symmetric two-body potential of 
short range, and g > is a coupling constant. The phase space of the coupled system is given by 

Ttot := ^particle X ^ BG (9) 

and its Hamilton functional by 

H(X, P; p, p) = H(X, P) + H BG (P, P) + ^(X; (3, P) (10) 
where if is as in ([T]), %bg as m © an d Hi as in ([8]). 



B. A list of models of decreasing difficulty 



(1) G-model : This is the model specified by ([9]) and (|10p . with p < oo, A > 0, g > 0. 



B A list of models of decreasing difficulty 
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(2) E- model : We consider the " Bogoliubov limit", 



In this limit. 



p — > oo, A — > 0, A/9 = const := - 
g — > 0, g-v/p =: ^, k and v fixed 



H(AT, P; (3, 0) = H(X, P) + H^tf, 0) + \X; (3, /3) 



(11) 



where 



n%l(Pj) = J d 3 x [^-|V/3(x)| 2 + K jd 3 y Stp{x) <&(* - y) »/3(y) 
(8, (3) = 2v J d 3 x W(X - x) 

It is useful to introduce real coordinates, 



(12) 



with 



= -j=(<p(x)+iir(x)) 



(13) 



Then 



«gg = J d 3 x [±- (Vvr) 2 (x) + (V^) 2 + d 3 y <p(x)*(x - y)<p(y)] (14) 



and 



Uf ] =V2vj d 3 x W(X - x)<p(a 



(15) 



The equation of motion for the Bose gas (g = 0) are 



1 



It follows that 



ip = {H,<p} = — Avr 
2m 

7T = VH,7r} = Aip + k$ * if 

2m 



(16) 



which is a wave equation for <p with constant coefficients that can be solved by Fourier transforma- 
tion. 



A Static solutions: X = P = = 
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One concludes that the frequency, £1, of a plane wave with wave vector k E M 3 is given by 



n «=i*iv(^) 2+ 



k y K$(k) 



2m/ 2m (17) 

Here, = i s the minimal propagation speed of a Goldstone mode (sound wave) in the 

Bose-Einstein condensate of the Bose gas. 

The equations of motion of the coupled system are given by 

MX = P, P= -VV{X) -u J d z x VW(X - x) 9fy3(x) 



ip( x ) = — — A/3(x) + - A ($ * K/3) (x) + vW x (x) 
2m 4 



(18) 



where W x {x) := W(x - X). 

(3) C-model : A = 0, p < oo, g > 0; see [4|. 

(4) B-model : A = 0, p — > oo, with g^fp = v fixed; see [fj, The equations of motion in 
the B-model are those given in (|18p . with k = 0. 



III. SOME SPECIAL SOLUTIONS OF THE EQUATIONS OF MOTION 
A. Static solutions: X = P = ft = 

Let X* £ IE 3 be a critical point of the external potential (with X* arbitrary if V = 0). Let /3* be a 
critical point of the energy functional 

£(X.;PJ) := J d 3 x [-L|V/3(x)| 2 + 5 W(X, -x)(|/3(x)| 2 + 2^^/3(x)) 

+ Xj d 3 y (|/3(x)| 2 + 2^-p K/3(x)) *{x - y) (|/3(y)| 2 + 2J~ P K/%)) 

A critical point /3* can be constructed with the help of variational calculus. Then 

(-) (-) \ 

X = X*,P = 0, /3 (x) = /3* (x) J is a static solution of the equations of motion. If X* is a non- 



degenerate local minimum of V one expects that there exist solutions of the equations of motion 
close to this static solution and converging to it, as time t tends to oo. 
In the E-model, /3* is the real solution of the elliptic equation 

+ ^-vW*-, (19) 



B Traveling waves for V = 
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which can be solved by Fourier transformation. If k<E>(0) = K J d?x <&(x) > then the solution 
exhibits exponential decay in x. 



B. Traveling waves for V = 

Ansatz 

X, = vt 4- Xn. Pt = Mv. v G R 3 ) 

(20) 



X t = vt + X , P t =Mv, v 6 
/3 t (x) = 7„(x - X t ) 



Let VF be smooth, spherically symmetric and of rapid decay at oo. We consider the E- model, with 

v * = V ^ai^ > ^' see Then, for \v\ < v c < the equations of motion (fTBj) have a traveling 

wave solution of the form (1201) . and 



7„ decays exponentially fast at oo, for \v\ < v c ; j v has power-law decay at oo, for \v\ = v c . 

As \v\ — > 0, 7„ converges to /3* (see (fT9|) ), with X* = Xo- If &(x) = 5(x) then v c = v*. 
The proofs of these claims follow from easy calculations, using Fourier-transformation. 
For \v\ > v c « v*, the equation 

k 2 k 

+ t*(*0 — « ■ A; = 

2m 4 v ; 

has real solutions, k = k v G M 3 , with fc„ || v. Then there do no£ exist any traveling wave solutions 
for V = 0. If the ansatz (I20p is plugged into the equations of motion one finds that P ^ 0, with 
P anti-parallel to v; i.e., a non-zero friction force, F v , arises. This force is caused by emission of 
Cerenkov radiation by the particle, and its precise form can be inferred from Fermi's Golden Rules. 
If K = (i.e., v* = 0) the expression for F v is very simple: 

F v = -const v-—^ ^ dp p\W(^p)\ 2 , (21) 
where the constant is proportional to v. The behaviour of \F V \ is indicated in the following figure. 



D Significance of special solutions 
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Our discussion suggests to proceed to study forced traveling waves, with V(X) = —F-X, for a 
non- vanishing external force FgK 3 . 



We study the motion of a tracer particle under the influence of a constant external force FgK 3 , 
with V(X) = —F ■ X. For simplicity, we consider the B- or the E- model. 

n 

Theorem 1. [S] Suppose that W is smooth, spherically symmetric and of rapid decay at infinity. 
Then there exists a positive constant F max < oo such that 

(1) if \F\ < F max the equations of motion U8\) of the B- or E-model have two traveling wave 
solutions of the form [2U\) with velocity v \\ F, \v\ > v c « v*; 

(2) if \F\ = F max there exists a unique traveling wave solution with v \\ F, \v\ > v c ; 

(3) if \F\ > F max there do not exist any traveling wave solutions of the form \2U\) . 

Remarks . To prove this theorem we must look for solutions of the equations of motion of the form 



If \v\ > v c the Fourier transform, %, of the second equation in (|22p is singular on a spherical surface 
in k-space, and it follows that the first equation in l\22\) does not have any solutions unless \F\ > 0, 
as discussed in section 1111 Bl 

The right side of the first equation in (|22p is a function of \v\ that vanishes, for \v\ < v c , tends to 
0, as \v\ — > oo, and has a unique maximum at some value, v* , of \v\. It is quadratic in W. 

From these observations the proof of the theorem stated above can easily be inferred. For \v\ > v c , 
\F\ ^ 0, the solution, 7„, of the second equation describes a conical wave whose tip is at the 
position, X t , of the particle. 



C. Forced traveling waves 



X t = vt + X , P t = Mv, Pt{x) = i v (x - X t ). 
The equations of motion (|18p then imply that 




(22) 



D. Significance of special solutions 



The construction of special solutions, as outlined in Sects IrTTAt IIII Bl and [ITTCl is elementary, and 
one may wonder why such solutions are of any interest, at all. The reason is that one expects 
that rather general classes of (time-dependent, non-stationary) solutions of the equations of motion 
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approach the special solutions constructed above, as time t tends to oo, (in a mathematically precise 
sense). We briefly discuss this expectation in several situations. 

(1) Let us suppose that V(X) has a non-degenerate (quadratic) minimum, V(X*) < 0, at a point 
X* £ E 3 , with V(X) — > 0, as \X\ — > oo. We consider the G-model and choose initial conditions, 
(Xq, Pq, fio), whose total energy, T~L(Xq, Pq; fio, (3q), is strictly negative. Then we expect that the 
solution, (Xt, Pt, fit), of the Hamiltonian equations of motion corresponding to the Hamilton func- 
tional H(P, X;fi,fi) in (fTUj) approaches a static solution, Xf — > X*, P t — > 0, fi t — )■ /3*, as t — )■ oo, 
where /?* is a minimizer of the energy functional £(X*; fi, fi) introduced in section llll A[ provided the 
coupling constant g is sufficiently small. (If the potential V is as specified above, but V(X) — > oo, as 
\X\ —7- oo, then one expects that the static solution (X*, P = 0, /?*) is asymptotically stable without 
further assumptions on initial conditions of finite total energy and on the coupling constant g.) 

(2) We consider the B- or E-model, with V{X) = 0. We choose initial conditions, (Xq, Pq, /3q), of 
energy T~L(Xq, Pq; /3q, fio) (see (I12j) ) not "much larger" than yi)^ (with v c = = 0, for the B-model) 
and g small enough. Then we expect that the solution, (Xt, Pt, Pt), of the equations of motion (|18D . 
with the given initial conditions, approaches a traveling wave solution (|20p . with \v\ < v c , as t — > oo; 
(convergence of fit is understood locally, in arbitrary balls of finite radius around the position, Xt, 
of the tracer particle) . 

(3) We consider the B- or E-model, but with V(X) = -F ■ X, < |F| < F max . Let u< > 
and Vy > v < be the two speeds for which the equations in (|22p have solutions, with |t;| = v < or 
\v\ = Vy. We first discuss the stability of solutions with |f| = v < . We expect that a solution 
of the equations of motion (|18p corresponding to an initial condition (Xq, Pq, fio), with Xq G E 3 
arbitrary, |Po| << M v <> %bg?(^0"^o) small enough, converges to a traveling wave solution of the 
type discussed in the theorem of section fill C\ with \v\ = v < , as time t — > oo, provided g is small 
enough. 

I ii contrast, if |Po| <i Mv>, Hbg(.Poi @o) sman > the solution (Xt, Pt, fit) of the equations of motion 
(|18p will in general not approach a traveling wave solution with \v\ = w>. Rather, we expect that 
\Pt\ diverges, as t — > oo. (Similar behaviour is expected if \F\ > F max .) 

Most of these expectations remain unproven, except for expectation (2) in the context of the B- 
model. 



IV. A THEOREM DESERVING THIS NAME 



In this section, we consider the B- and C-models, which are the only models for which non-trivial 
results have been established, so far; see ((J, 0|,0|. We state our main result for the B-model; but 
a very similar result has also been proven for the C-model. We propose to study solutions of the 
equations of motion (I18p . for k = 0. 

We choose units such that M = m = 1. The two-body potential W is assumed to be smooth, 
spherically symmetric and of rapid decay at oo, with |IT^(0)| = (2-7r) _3 / 2 |/ d?x W(x)\ = 1. The 
wave field fi is normalized appropriately. The external potential V is assumed to vanish identically. 
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Then the equations of motion (|18p take the form 

X t = Pt, Pt = v J d 3 x (VW Xt )(x) m(x) 



»ftO»0 = ~\{Wt){x) + W x *(x) 



(23) 



where W Xt (x) := W(x — X t ), v = 0(1) a constant. Note that the equations (f2"3"j) have static 
solutions 

X t = X*, P t = 0, k(x) = /3,(x) = 2{&- l W x *){x) 

indexed by X*, and we propose to analyze their asymptotic stability. Our main result is the 
following theorem. 

Theorem 2. |J/ Under the assumptions described above, there exists an interval I D [0, 0.66) D 
[0, ^] swc/i i/tai, /or any 5 £ I, there is some 6q = eo(5) > with the property that if 



Pol \\(x) 3 Po\\2, ||V/3 || 2 <e 



where (x) := ^Jl + |x| 2 , t/ien 



/or some finite constant c, and 



Pt\ < c(t)~5~ S , as t -> 00 (24) 



lim \\/3 t - 2A- 1 W Xt \\ oa = 

i— >oo 



Remarks: 



[1) If 5 is chosen to be larger than ^ (note that I contains such values of 5) then 

Xt = Xq + P s ds — > X*, as t — > 00 
Jo 



with < 00. 

(2) We suspect that \Pt\ ~ t~z~ & \ as i — )• 00, for an exponent 5 > | independent of initial 
conditions, as long as |Po|j IK^^Ajlk) ||V/3o||2) are finite. However, this has not been proven, so 
far. (If the condition that ||(x) 3 /3o||2 < 00 is dropped then 5 may depend on initial conditions. 
For, consider the equation 

v t = F vt (25) 

with F v as in (|2ip . This corresponds to /3t(x) = j Vt (x — Xt). The solution of (|25p behaves like 
\vt\ ~ (t) -1 , i.e., 5 = 2- However, ||/3t||2 is divergent!) 



(3) The theorem stated above has interesting applications in the analysis of decoherence in 
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systems of a heavy quantum-mechanical particle moving through a Bose-Einstein condensate. 

Although the initial moves in the proof of the theorem described above are quite natural and 
fairly standard, the technical details are suprisingly tedious and involve exhibiting very subtle 
cancellations among terms that do not, a priori, exhibit adequate decay, as time t — > oo, in order 
to "close the estimates". 

It would be of considerable interest to study systems of finitely many tracer particles moving 
through a Bose-Einstein condensate, or a gas of such particles. It is not very hard to come up with 
plenty of interesting and plausible looking conjectures concerning the behaviour of such systems. 
But it appears to be very difficult to rigorously prove these conjectures. 

Note : For further results on friction, see also [H, 0, 0| • 
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